VARIABILITY OF SOUTHERN AND NORTHERN 

SKR PERIODICITIES 



L. Lamy* 



O 

^ Abstract 

I^T j Among the persistent questions raised by the existence of a rotational modu- 

lation of the Saturn Kilometric Radiation (SKR), the origin of the variability of 
the 10.8 hours SKR period at a 1% level over weeks to years remains intriguing. 
While its short-term fluctuations (20-30 days) have been related to the variations of 
p 1 the solar wind speed, its long-term fluctuations (months to years) were proposed to 

PQ be triggered by Enceladus mass-loading and/or seasonal variations. This situation 

^ has become even more complicated since the recent identification of two separated 

Qh periods at 10. 8h and 10. 6h, each varying with time, corresponding to SKR sources 

Q located in the southern (S) and the northern (N) hemispheres, respectively. Here, 

^ six years of Cassini continuous radio measurements are investigated, from 2004 (pre- 

^ equinox) to the end of 2010 (post-equinox). From S and N SKR, radio periods and 

"—I phase systems are derived separately for each hemisphere and fluctuations of radio 

^ periods are investigated at time scales of years to a few months. Then, the S phase 

>► is used to demonstrate that the S SKR rotational modulation is consistent with an 

intrinsically rotating phenomenon, in contrast with the early Voyager picture. 

o 
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^ 1 Introduction 



The Saturn Kilometric Radiation (SKR) is an intense non-thermal radio emission pro- 
duced by auroral electrons moving along magnetic field lines, predominantly on the dawn 
sector [Kaiser et al., 1980]. Its regular pulsation, whose origin still remains unexplained, 
was originally interpreted as a clock-like rotational modulation triggered by the planetary 
magnetic field, and thus directly relating to the planetary interior. The Voyager deter- 
mination of the SKR period (10h39min24±7s or 10.657±0.002h) [Desch and Kaiser,1981] 
was adopted as the official rotation rate of the planet by the International Astronomical 
Union [Seidelmann et al., 2001], and used to define the Saturn Longitude System (or SLS, 
hereafter named SLS 1). In this system, the SKR occurrence or intensity is organized by 
sub-solar longitudes and peaks at a fixed arbitrary value (historically set at 100°). 
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Ulysses observations later showed that the observed SKR period (^10. 8h) is not con- 
stant but varies with time by 1% over years [Galopeau and Lecacheux, 2000], which 
was confirmed by Cassini measurements [Gurnett et al., 2005]. Further studies identified 
short-term 20-30 days oscillations of the southern SKR period [Zarka et al., 2007] corre- 
lated to the fiuctuations of the solar wind speed (while the SKR intensity is correlated to 
the solar wind dynamic pressure) supporting previous modeling work [Cecconi and Zarka, 
2005], together with long-term yearly oscillations also found in magnetospheric plasma 
and magnetic field data and attributed to Enceladus mass-loading [Gurnett et al., 2007]. 
New longitude systems (SLS 2,3), based on polynomial fits of the SKR long-term phase 
variations successively replaced SLS 1 [Kurth et al., 2007, 2008]. 

The identification of a second SKR period around 10. 6h [Kurth et al., 2008] has strong 
implications on the validity of the longitude systems described above. Indeed, radio 
periods at ^10.8h and ^10.6h were respectively attributed to SKR emissions emanating 
from the southern (S) and the northern (N) hemispheres [Gurnett et al., 2009]. It is 
interesting to note here that Voyager essentially observed northern SKR, while its orbital 
motion out of the ecliptic enabled Ulysses to observe both hemispheres close to Saturn's 
equinox. As both Cassini-derived SKR periods vary with time and come close together 
about Saturn's equinox (Aug. 2009), these authors proposed that seasonal variations of 
the solar illumination may induce a different magnetosphere-ionosphere coupling in each 
hemisphere, and trigger different S,N radio periods varying in opposite ways, with an 
expected crossing about equinox, when the sun illuminates both hemispheres equally. 

Section [2] describes the employed radio dataset and the periodogram technique used to 
derive separate S and N SKR periodicities as a function of time. Section [3] investigates 
fiuctuations of these periods at time scales of years to a few months, and implications on 
magnetospheric dynamics. Long-term phase systems are then built for each hemisphere. 
Finally, section [4] investigates the organization of S SKR sources as a function of S phase. 



2 Radio observations and harmonic analysis 
2.1 Northern and southern SKR emissions 

Taking advantage of quasi-continuous SKR observations acquired by the Cassini Radio 
and Plasma Wave Science (RPWS) experiment since 2004, radio data were processed as 
detailed by Lamy et al. [2008a] to obtain regular dynamic spectra of circular polarization 
degree V and radiated power P (in W.sr~^) from 3.5 to 1500 kHz between 1 January 2004 
(DOY 1) and 19 October 2010 (DOY 2484) with a 3 min time resolution. A time series 
of the total radiated power Pskr{^) was then obtained by integrating P over the spectral 
range 40 to 500 kHz, that corresponds to well-defined SKR bursts and avoids narrowband 
low frequency emissions below 40 kHz, that display periods comparable to SKR ones but 
with a significant phase shift [Ye et al., 2010]. 

Southern (S) and northern (N) radio emissions were separated on the basis of (i) the SKR 
property to be primarily emitted on the extraordinary mode [Kaiser et al., 1984; Lamy et 
al., 2010 and references therein], that displays left-handed (LH, V > 0) and right-handed 
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(RH, V < 0) circular polarization for S and N sources, and (ii) its visibility domain, 
that mainly illuminates its hemisphere of origin and extends down to 20° latitude in the 
other hemisphere (see Fig. 11 of Lamy et al., [2008a]). S and N SKR emissions were 
thus identified by LH and RH emissions respectively observed from latitudes A^c < 20° 
and \sc > —20°, then integrated from 40 to 500 kHz to build time series ^skr,s{^) and 
^ skr,n{^)' This technique permits one to consider all near-equatorial observations where 
S and N SKR are observed together (most of Cassini measurements), generally excluded 
by separating S,N SKR from a single geometrical selection on the spacecraft latitude. 
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Figure 1: Lomb-Scargle 2D periodograms of (a) log PsKR(t), (b) log PsKR,s(t) ^.nd (c) 
log PsKR,N(t) SiS a function of period and time. Individual power spectra (columns, in arbi- 
trary units) were computed everay day over a 200-days long window. Solid yellow lines show S 
and N periods, as derived from panels (b) and (c). 
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2.2 Lomb-Scargle normalized periodogram 
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The Fast Fourier Transform (FFT) is the technique employed most often to perform the 
spectral analysis of regular time series. However, it cannot deal with unevenly spaced 
data and, when applied to discontinuous series re-interpolated on a regular basis, it yields 
noisy power spectra with a significant spectral leakage. The Lomb-Scargle (LS) normalized 
periodogram is a powerful alternate technique directly applicable to irregularly sampled 
data that possesses useful properties [Lomb, 1976; Scargle, 1982; Horne and Baliunas, 
1986]. The periodogram analysis is equivalent to least-squares fitting of weighted sine 
curves to the data. Its normalization by the total variance leads to less noisy power spectra 
and, while the FFT spectral resolution is fixed by the time step of the original time series, 
the periodogram can be oversampled in order to improve the determination of significant 
spectral peaks. Due to non strictly continuous RPWS measurements together with the 
applied data selection, regular SKR power series described above intrinsically contain 
unevenly spaced non- null signals, for the spectral analysis of which the LS technique is 
particularly adapted. 

3 Variability of SKR periodicities 
3.1 Long-term variations 

Figure [l] displays 2D periodograms of the LS power spectra of log Pskr{^)^ log Pskr,s{^) 
and log Pskr,n{^) as a function of period and time, computed with a 200-days long sliding 
window and a time resolution of 1 day. Figure [l^ reveals the two main SKR periods around 
10. 8h and 10. 6h, each slowly varying with time, corresponding to southern and northern 
SKR sources. This is illustrated in Figures [1)3 and[T]3, where S and N SKR mainly pulse at 
periods about 10. 8h and 10. 6h, identified by solid yellow lines hereafter labelled T^,A^(t). 
In Figure [T^, the intensity of the S peak is dominant over 2004-2007 and then decreases 
relative to the intensity of the N peak when approaching the equinox. This refiects the 
dominant observed signal which corresponds, for comparable observing conditions of both 
hemispheres (for instance at the equator), to the intrinsic SKR power radiated by each 
of them, itself observed to change with seasons: while N SKR was more intense during 
northern summer at the Voyager epoch [Kaiser et al., 1984], S SKR was predominant 
during the first years of southern summer observed by Cassini [Lamy et al., 2008a]. 

Ts{t) and T]sf{t) are plotted together in Figure [2] (black lines). They display clear opposite 
trends, with a correlation coefficient reaching c = —0.95, maximal for a lag of days. Both 
periods cross about 7 April 2010 (DOY 2050), namely approximately 8 months after the 
equinox. The S period is maximal around 16 June 2007 (DOY 1263) and the N one after 
2006 is minimal around 20 October 2007 (DOY 1389). Out of this interval, both periods 
significantly shift from their extrema. 

This anti-symmetrical behaviour supports a seasonal origin. However, the delay of 8 
months between the equinox and the actual reversal of both periods indicates a non-linear 
response of the magnet osphere. Focussing on seasonal variations of the solar illumination. 
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Figure 2: S and N SKR periods derived from a Lomb-Scargle analysis of SKR power using 
200-days long (black, Figures^,c) and 100-days long (blue) sliding windows. The equivalent 
resolution of FFT is indicated by double arrows on the left side of the Figure. Both periods 
cross at day 2010-097 (DOY 2289). Orange dashed curves indicate the sub-solar latitude Xgun 
and -Xsun (right axis) that cross at the equinox of 11 Aug. 2009 (2009-223 or DOY 2050). 

Gurnett et al. [2009] proposed that the temporal variation of Saturn's inchnation controls 
the atmospheric Pedersen conductivity in each hemisphere, which in turn triggers the 
strength of associated auroral field-aligned current systems, and thus the intensity of the 
torque (opposite to corotation) at the footprint of the magnetic field lines where these 
currents close. As the Pedersen conductivity varies quasi-instantaneously with the solar 
illumination, this model needs to be augmented with a significant mechanical slippage 
responsible for the observed lag. Also, while they are close to each other (within 100 days), 
the extrema of both periods are shifted with respect to the ones of the solar illumination 
(dashed lines in Figure [2]) by approximately 2 years, suggesting that additional effects 
shall aflFect the variation of radio periods. 

Interestingly, N SKR displays a secondary peak at the S period (obvious in 2005 and 2007). 
Although it cannot be formally excluded, it is unlikely that such an intense secondary 
peak results from a non-ideal data selection with a possible residual contribution of S SKR 
sources to Pskr,n{^)- Moreover, this secondary peak remains when separating S and N 
emissions from observing latitudes above thresholds of 10° or 20°, and is supported by 
equatorial observations (see example in 2007 displayed by Figure ^ , where part of the N 
emission pulses together with the S one. This result is thus likely to refiect a real physical 
dual modulation. This is not surprising since SKR sources have been identified on closed 
field lines [Lamy et al., 2010; Bunce et al., 2010], and that auroral electrons accelerated 
in one hemisphere could ultimately reach the other one, as proposed for interpreting Jo's 
multiple footprints [Bonfond et al., 2008]. Moreover, both periodicities have been detected 
in northern electron observations [Carbary et al., 2009], as well as in equatorial magnetic 
field oscillations [Provan et al., 2011]. 

In a parallel paper, Gurnett et al. [2010] specifically investigate the crossing of SKR 
periods. They show that periodicities observed by Ulysses (determined by FFT) are 
consistent with the observation of S,N hemispheres, and crossed about 9 months after the 
equinox of 1995. With a diflFerent technique of period determination, they also identified 
a crossing of both SKR periods from Cassini data approximately 7 months after the 



6 L. Lamy 

equinox of 2009. Although they are derived with a 240-days shding window excluding 
near-equatorial observations, these periods are roughly consistent with the ones derived 
above. However, the 8 months delay measured from the present LS analysis is in better 
agreement with the one derived from Ulysses data. Also, the final S SKR period given by 
combining Ulysses/Cassini data displays significant long-term changes shifting from the 
simple smooth evolution of the planet inclination. 
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Figure 3: (a) Plot of log Pskr,s (^) normalized to the average over one rotation as a function of an 
arbitrary phase computed from a fixed reference period of 10. 7928 h (enabling direct comparisons 
with Kurth et al. [2008]), between 2004-001 and 2010-193. Each rotation is replicated four times 
along the y axis for clarity. The S SKR rotational modulation (given by^s(t) — 0° j is displayed 
in solid yellow, while the SLS 3 one (given by sub-solar longitudes = 100°^ is shown in solid light 
blue. Both correctly track S SKR maxima, (b) Identical plot for log Pskr,n (t) ^nd an arbitrary 
reference period of 10.6 h. (c) Cassini orbital parameters (Xgc the latitude, LTgc the local time 
and Tsc the distance to the planet ^s center in kronian radii Rs, with 1 Rs — 60268 km). 
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3.2 Phase systems 

Using continuous SKR periods, a separate phase system was built for each hemisphere 
from 1 January 2004 (DOY 1) to 12 July 2010 (DOY 2385). The phase of the SKR 
rotational modulation $5',Ar(t) was numerically integrated from: 

/360 
Tf^^^dt + ^o^s^N (1) 

where t is the time and ^o,s,n an arbitrary reference. As $5',Ar(t) is defined modulo 
360°, SKR maxima occur at a fixed phase. ^o,s,n was chosen to reference SKR peaks at 
^s,N{t) = 0°, by fitting log Pskr,s{^s) and log Pskr,n{^n) (respectively averaged over 
the intervals 2004-2008 and 2005-2008, where rotational modulations are well defined) 
with cosine functions. 

Figure [3] illustrates the temporal evolution of S,N SKR phase drifts with respect to ar- 
bitrary constant rotation periods. Yellow solid lines show that the above defined SKR 
phases correctly track long-term SKR maxima for each hemisphere, with more noisy in- 
tervals after the equinox or at the end of 2007. However, the latter intervals correspond 
to near-equatorial observations from the dusk sector, from which the visibility of dawn 
radio sources can be questioned (see section ^ , as suggested by previous modeling work 
[Lamy et al., 2008b] and direct observations [Cecconi et al., 2009]. ^^(t) compares to the 
sub-solar longitude defined in previous SLS models. In Figure $^(t) = 0° and SLS 
3 sub-solar longitude = 100° (solid light blue) match on the interval where the latter is 
defined, and differ by generally less than 20° to 30°. The relevance of the present phase 
systems is illustrated diflFerently in Figure [4] for 3 days of equatorial observations over 
which S and N SKR bursts (well identified in Figure |4j3) match $^,A^(t) = 0° (arrows). 

The S and N phase systems have been used to organize successfully the power radiated by 
UV aurorae [Nichols et al., 2010] in both hemispheres, which demonstrated the existence 
of a SKR-like diurnal modulation of atmospheric aurorae, which had been unsuccessfully 
searched for a long time. In parallel, Andrews et al. [2010] identified two systems of 
magnetic field oscillations, interpreted as independent high-latitude field-aligned current 
systems, that pulse at S and N SKR periods to within 0.01%, with upward current layers 
matching SKR intense dawn sources during S and N SKR peaks. In addition, Provan et 
al. [2011] showed that part of the 'jitter' observed in equatorial magnetic oscillations can 
be explained by beating effects between both periods. 

3.3 Mid- term variations 

To investigate fiuctuations of SKR periods at shorter time scales, the above LS analysis 
was applied on the same dataset with a 100-days long sliding window. The resulting 
S,N periods are plotted in blue in Figure [2] Although less continuously defined, they 
display clear fiuctuations of the order of a few months. Their correlation coefficient 
reaches c = —0.41, which increases to c = +0.27 once long-term trends (black) are 
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subtracted. Considering the limited number of oscillations, their approximately one-to- 
one correspondence supports an overall correlation. 

A symmetrical behaviour of T^,A^(t) may result from either orbital effects affecting the 
visibility of S and N SKR sources (whose intensity displays a clear LT dependence) or 
from a physical process affecting both hemispheres similarly. Orbital effects of the order of 
15 to 30 days clearly appear in Figures |3^,b and prevent to track 20-30 days fluctuations 
previously identified by Zarka et al. [2007] in pre-SOI observations. However, no obvious 
orbital variations of the order of a few months are visible, which makes plausible a similar 
control of both SKR periods by a common physical cause. Checking the possible presence 
of similar fluctuations in other magnetospheric periodicities will help to determine their 
origin. 




Figure 4: Dynamic spectra of (a) Eux density (normalized to lAU) and (b) degree of circular 
polarization for equatorial observations in 2007. LH (black) and RH (white) polarized bursts 
show the diurnal modulation of S and N X mode SKR. Top black arrows indicate SLS 3 SKR 
peaks. Bottom black and white arrows indicate S and N SKR peaks given by ^s,N(t) — 0°. 



4 Nature of the SKR diurnal modulation 

While most of periodic features observed in the kronian magnetosphere correspond to 
rotating phenomena, the SKR rotational modulation has been interpreted since the Voy- 
ager epoch as a strobe-like feature pulsing on the dawn sector. However, the radio source 
localization technique developed from the goniopolarimetric capability of Cassini radio 
instrumentation [Cecconi et al., 2009] has provided a new way to check this assumption. 
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From this method, Lamy et al. [2009] have shown that SKR sources actuaUy exist at aU 
LT, forming a statistical radio oval conjugate to the UV one and whose intensity strongly 
varies as a function of LT, peaking at 08:00 LT. This result challenges the Voyager picture. 

Here, ^^'(t) was used to organize the locus of S SKR sources (the most intense ones) 
derived from 6 years of individual 3- antenna measurements. They are less numerous than 
2-antenna measurements, but yield unambiguous wave directions, and therefore accurate 
source locus, whatever the wave polarization (for details, see Fischer et al. [2009]). The 
spectral range was limited to the 100-400 kHz SKR spectral peak, for which intensities at 
different frequencies are both high and comparable. Figures [5^ and [5)3 respectively show 
the median intensity and the occurrence of S SKR sources as a a function of their LT and 
^s- Figure [5]3 displays the associated orbital coverage with the occurrence of SKR sources 
organized as a function of the LT of the spacecraft and Uncertainties on the radio 
sources position have not been taken into account in Figures [5^,b, that merely displays 
an average behaviour. 

Figure [5]3 indicates that the S hemisphere was mainly visible when Cassini was between 
23:00 and 11:00 LT, while Figure [5^ illustrates that SKR sources were detected at all lon- 
gitudes, with a clear LT dependence of their intensity, peaking around 08:00 LT. Then, 
whereas a clock-like pulsation in Figure [5^ would correspond to highest intensities con- 
centrated around $5' = 0°, SKR sources surprisingly display a clear linear organization 
with LT($^), not seen in the orbital coverage (Figure [5]3), and consistent with an intrin- 
sically rotating phenomenon, materialized by the dashed arrows. But although the lowest 
intensities peak around $^ = 180°, as expected, the highest ones are observed around 
$5 = 300°, instead of 0°. This shift must be analyzed carefully as the intensity of de- 
tected sources merge contributions of a rotating feature, a strong LT dependence as well 
as important external effects (solar wind dynamic pressure, plasmoid activity) that can 
modulate the SKR flux by orders of magnitude. A complementary information is given 
by the occurrence of S SKR sources in Figure [5)3, less sensitive to overall modulations, 
which peaks well around = 0°. The low statistics for dusk SKR sources is likely due 
to the poor dusk orbital coverage, together with the result of the point-source assumption 
used in the goniopolarimetric analysis (if a dawn bright source and a dusk faint source 
are observed together, the derived direction will correspond to the most intense - dawn 
- source). Finally, several linear features are present both left and right to the dashed 
arrows, suggesting that the longitudinal extent of the rotating pattern is significant and 
covers more than one quadrant in LT. 

While this result needs to be separately confirmed by the symmetrical organization of 
SKR N sources (for which much less data fit the selection criteria before the equinox, both 
because of fainter intensity and poor visibility conditions) and the statistical analysis of 
2-antenna measurements, it departs from the Voyager picture and reconciles the dynamics 
of SKR sources with the ones of other periodic magnetospheric phenomena. The above 
phase systems are reliable until the equinox in the sense that $6',Ar(t) = 0° indicate the 
transit of S,N rotating features in the dawn sector, about 08:00 LT. Indeed, the good 
organization of SKR maxima observed in Figure |3^,b before mid-2009 suggests that, 
despite severely changing visibility conditions, auroral radio sources of the active morning 
region have been detected by Cassini for most of the mission. 
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Figure 5: (a) Median flux density (normalized to lAU) and (b) occurrence of S SKR sources 
organized as a function of their LT and (1^ bins of 15 min by 5^) from 2004-181 to 2010- 
193. The source LT was obtained from the radio source localization technique [Cecconi et al, 
2009] applied to 3-antenna measurements for Vsc < 201?5' and frequencies between 100 and 
400 kHz. This dataset corresponds to 231 observing days, to which has been applied a standard 
data selection [Lamy et al, 2010] to sources with footprint latitudes between —85° and —60° 
(surrounding the S radio auroral oval). These selection criteria were satisfied by 1% of the data, 
while the contribution of S sources detected after the equinox was negligible. The black dashed 
arrows materialize a source rotating at the S period, with a maximal intensity at 08:00 LT 
reached at = 0°. (c) Similar to (b) but as a function of the LT of the spacecraft and 

5 Conclusions 

In this article, accurate SKR periods were derived for each hemisphere from a 6-years 
long Cassini radio dataset analyzed by the Lomb-Scargle periodogram technique. S and N 
periods display clear anti-symmetrical yearly variations, consistent with a seasonal effect, 
together with symmetrical shorter term oscillations of the order a few months, whose origin 
remains to be identified. The continuous determination of S and N long-term varying 
periods, crossing 8 months after the equinox, was used to build phase systems for both 
hemispheres, vahd from 1 January, 2004 (DOY 1) to 12 July, 201 (DOY 2385), and freely 
available on the Cassini/RPWS/HFR data server: |http://www.lesia.obspm.fr/kronos| 

The S phase has been used to organize the locus of S SKR sources, which reveals an 
intrinsically rotating phenomenon. However, contrary to Jupiter, the search-light rota- 
tional modulation is aflFected by a strong LT dependence, which explains why the clock-like 
picture has remained a good approximation of the SKR periodic behaviour since Voyager. 
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